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A large body of evidence from lattice calculations indicates that instantons play a major role in the physics 
of light hadrons. This evidence is summarized, and recent results concerning the instanton content of the SU(3) 
vacuum, instanton contributions to the static potential, and a new class of instanton solutions at finite temperature 
are reviewed. 



1. INTRODUCTION AND PHYSICAL 
PICTURE 

1.1. Physics motivation 

One of the great advantages of the path in- 
tegral formulation of quantum mechanics and 
field theory is the possibility of identifying non- 
perturbatively the stationary configurations that 
dominate the action and thereby identify and un- 
derstand the essential physics of complex systems 
with many degrees of freedom. Thus, the discov- 
ery of instantons in 1975 |[) gave rise to great 
excitement and optimism that they were the key 
to understanding QCD. Indeed, in contrast to 
other many body systems in which the quanta 
exchanged between interacting fermions can be 
subsumed into a potential, it appeared that QCD 
was fundamentally different, with topological ex- 
citations of the gluon field dominating the physics 
and being responsible for a host of novel and im- 
portant effects including the 6 vacuum, the ax- 
ial anomaly, fermion zero modes, the mass of the 
77', and the chiral condensate. However, despite 
nearly a quarter of a century of theoretical ef- 
fort, it has not been possible to proceed analyt- 
ically beyond the dilute instanton gas approxi- 
mation In the intervening years, the instan- 
ton liquid model ||-|| provided a successful phe- 
nomenology and qualitative physical understand- 
ing, but a quantitative exploration of the role of 
instantons in nonperturbative QCD has had to 
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wait until lattice QCD became sufficiently sophis- 
ticated and sufficient resources could be devoted 
to the study of instanton physics. 

In recent years, significant progress has been 
made by a number of groups in understanding 
instanton physics on the lattice, and this review 
summarizes some of the highlights. This field 
is both exciting and frustrating. It is exciting 
because, at last, a substantial body of evidence 
shows that the original vision of the role of in- 
stantons in QCD is borne out in nature. How- 
ever, as will be clear below, it is also frustrating 
because of intrinsic limits to the precision with 
which one can apply ultimately semiclassical con- 
cepts. Since there is more interesting material 
than could be accommodated, I will omit signifi- 
cant developments in spectral methods, which are 
covered in another talk and recent work on 
abelian projection and monopoles, a subject to 
which I am not able to provide any additional 
insight. 

1.2. Aspects of continuum instanton 
physics 

To put lattice investigations in context, it 
is useful to recall relevant aspects of contin- 
uum instanton physics. Working in Euclidean 
time, we evaluate a path integral of the form 

Hence, as in statistical me- 
chanics, the weight of a configuration depends not 
only on its energy, but also on its entropy - the 
number of ways it can be realized. In addition, 
tunneling solutions arise as periodic classical so- 
lutions in an inverted potential. Vortices in two- 
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dimensional spin systems provide a useful anal- 
ogy. Classified by winding number, the cost of the 
localized core energy is compensated by entropy 
so that vortices can become the dominate degrees 
of freedom as in the Kosterlitz-Thouless phase 
transition. Similarly, we may expect BPST in- 
stantons Q which connect degenerate minima of 
differing winding number with the self-dual gauge 
potential A a ^{x) — -j^ffi having scale-invariant 

action S = \J d i xF^F£ l/ = %-K 2 g- 2 = S to 
have a significant presence in the vacuum due to 
the high entropy associated with translation, size, 
and color orientation. The action and topological 
charge density are localized around the center as 
± V£ = = tM^Y ■ T J e tunneling 

rate is dm ~ (^-) 2N ^d 4 x ( ^e~^^ (Ap) 11 ^ , 
where the prefactor and the running of the cou- 
pling constant in the last factor produce a dis- 
tribution of instantons ~ p 6 for SU(3). Physi- 
cally, we expect this distribution to be cut off at 
large p by interactions between instantons and by 
fluctuations when the amplitude of a sufficiently 
large instanton becomes small relative to quan- 
tum fluctuations. It is the difficulty in treating 
these infrared effects that has stymied analytic 
progress. 

From the axial anomaly, 
d» £/ ^7m7sV> = 2m ^ + ^F^F^, 
the topological charge satisfies the index theorem 
and, for periodic systems, may be expressed in 
terms of fermion eigenfunctions 

Q=^fFF = n L -n R = mEx I ^:t i :\ 
where til and tir denote the number of fermion 

zero modes. For an isolated instanton, the zero 

mode is ij} a (x) = afff+^p^ 0- In the limit of 

light quarks, the Greens function for Nf quarks 

reduces to the product of zero modes 

and gives rise to the 't Hooft interaction. Thus, 
light quarks propagate by zero modes which in 
turn arise from instantons. Based on large N ar- 
guments, the Veneziano-Witten formula |?],|| re- 
lates the rj mass to the topological susceptibility 
in the pure gluon sector 

X = J ^(Q(x)Q(0)) = Sj(m* + m%, - 2 m \) 
yielding the expectation that \ = (180 Me V) 4 . 



Finally, for an ensemble of instantons and anti- 
instantons, we expect mixing of their zero modes 
to generate a finite density of states near zero vir- 
tually. Hence, by the Banks-Casher relation jo), 
-(H>) = JdXp(X)j^ —,^ = 0), they ac- 
count for the chiral condensate. Subsequent sec- 
tions will discuss the extent to which this contin- 
uum instanton physics can now be seen on the 
lattice. 

1.3. The instanton liquid 

The instanton liquid model [^|-|| provides an 
economical phenomenology of instanton mediated 
quark propagation in the QCD vacuum. The 
integral over all gluon fields that one evaluates 
in lattice QCD using an ensemble of configura- 
tions sampling the action is replaced by an en- 
semble of instanton and anti-instanton configura- 
tions. Although the model has been refined to 
account for interactions between instantons and 
the fermion determinant, the gross features can 
be seen for the case of an ensemble of instan- 
tons and anti-instantons of size p ~ 1/3 fm and 
density n ~ 1 fm 4 randomly distributed in space 
and color orientation, where the values of p and n 
are determined from the physical gluon and chiral 
condensates. 

Due to the opposite shifts of left and right- 
handed fermion levels in the presence of an in- 
stanton, one may think of the 't Hooft interac- 
tion as a vertex that absorbs left-handed par- 
ticles of each flavor and creates corresponding 
right-handed particles, and vice versa for anti- 
instantons. Mesons then propagate in the QCD 
vacuum by the hopping of quark-antiquark pairs 
between these vertices, and the qualitative fea- 
tures of the channel dependence arises naturally. 
Considering two flavors for simplicity, a spin-zero 
pion propagates by a u^dL pair interacting with 
an anti-instanton to produce a u l^r pair which 
in turn interacts with an instanton to return to 
a Uftdh pair. Since the 't Hooft vertex is most 
attractive in this channel and the interaction can 
act in all orders, the pion is the most strongly 
attractive meson channel. In contrast, the spin- 
one rho propagates only in second order, with an 
anti-instanton taking a u udji pair to a u^dLdudii 
state which is then returned to a upjdn pair by 
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an instanton, so we expect the interaction in this 
channel to be much weaker. Similarly, the scalar 
meson channel is repulsive and one expects the 
nucleon channel to be more attractive than the 
delta channel. The chiral condensate arises nat- 
urally in this picture by the fact that the zero 
modes for isolated instantons mix in the instan- 
ton liquid giving rise to a finite density of states 
at low virtuality. 

The tendency of instanton-anti-instanton pairs 
connected by light quarks and antiquarks to form 
dipoles oriented in the thermal direction sug- 
gests a possible mechanism for chiral symmetry 
restoration at high temperature and may be man- 
ifested by dipoles wrapped around the periodic 
thermal direction above Tq. 

2. EVIDENCE FOR INSTANTONS IN 
HADRON STRUCTURE 

The expectations raised by the our analytical 
understanding of the role of instantons in QCD 
and the physical picture of quarks hopping be- 
tween zero modes associated with instantons have 
now been borne out by a broad range of lattice 
calculations as described below. 

2.1. Vacuum current correlation functions 

Vacuum correlation functions for space-like 
separated hadron currents calculated in lattice 
QCD display the qualitative behavior expected 
from the 't Hooft interaction and agree semi- 
quantitatively with the instanton liquid model. 
As emphasized in ref |icj , correlation functions 
of the form R(x) = (0|TJ M (x) J M (0)|0> charac- 
terize the spatial and channel dependence of the 
interaction between quarks and antiquarks and 
thus supplement hadron bound state properties 
like phase shifts supplement deuteron properties 
in characterizing the nuclear interaction. The ra- 
tio R(x) / Rq(x) of the interacting to free corre- 
lator has been calculated in quenched QCD for 
the following meson and baryon currents p| , 
J = u'j^d, U7 M 7 5 c?, U75<f, ud, e a bc[c a Cj tl u b ]'y tl ~/5d c 1 
and e a bc[u a C'yfj,u b }u c . Results are consistent with 
dispersion analysis of c + -c~ and other data in rel- 
evant channels, and typical results are shown in 
Fig. 1. 
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Figure 1. Vector (V) and Pseudoscalar (P) cor- 
relation functions from lattice calculations 111] 
(solid circles), from the random instanton model 
p^ , p^| (open circles) and from dispersion analy- 
sis p0| (long dashes). The solid curve fit to the 
lattice data is made up of a continuum contribu- 
tion (short dashes) plus a resonance contribution 
(dotted curve). 

Note for subsequent reference, that the solid 
curve fit to the lattice data may be decomposed 
into a continuum contribution concentrated near 
the origin (short dashed line) and a resonance 
contribution arising from the rho or pion which 
dominates in the region of 1 fm and beyond (dot- 
ted line). Quenched calculations at /3 = 6.2 |14| 
corroborate the original (3 — 5.7 results. 

2.2. Comparison of results with all gluons 
and only instantons 

One dramatic indication of the role of in- 
stantons in light hadrons is to compare observ- 
ables calculated using all gluon contributions 
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Figure 2. Vacuum correlator in the rho channel 
calculated with all gluons (upper left) and with 
only instantons (upper right) ||l7f . The rho meson 
ground state density-density correlation function 
calculated with all gluons (solid circles) and with 
only instantons (open circles) is shown below. 



wave functions of charmed and B mesons would 
be drastically changed. 

As shown in Fig. ^, however, the properties 
of the rho meson are virtually unchanged. The 
vacuum correlation function in the rho (vec- 
tor) channel and the spatial distribution of the 
quarks in the rho ground state, given by the 
ground state density-density correlation func- 
tion Jl8| (p\qjoq(x)qjoq(0)\p), are statistically in- 
distinguishable before and after cooling. Also, as 
shown in Ref. [l7| , the rho mass is unchanged 
to within its 10% statistical error. In addition, 
the pseudoscalar, nucleon, and delta vacuum cor- 
relation functions and nucleon and pion density- 
density correlation functions are also qualitatively 
unchanged after cooling, except for the removal 
of the small Coulomb induced cusp at the ori- 
gin of the pion. Similarly, the axial charge ma- 
trix elements specifying the spin content of the 
nucleon, (PS\(qYij 5 q)\PS) = 2S'Aq, are quite 
similar when calculated with all gluons and only 
instantons fll9|] . This result suggests that when 
all the quark and gluon contributions are calcu- 
lated, instantons, which are the natural mecha- 
nism to remove helicity from the valence quarks 
and transfer it to gluons and sea quarks, may ac- 
count for the proton "spin crisis" . 



with those obtained using only the instantons re- 
maining after cooling. The basic idea of cool- 
ing jl5];[l6| i an d the variants described in Section 
3, is to iteratively locally relax a lattice gluon con- 
figuration and thereby approach the nearest local 
minimum of the classical action. Thus, one may 
extract the instanton content by removing local 
fluctuations while preserving topologically stable 
excitations. When a (3 = 5.7 quenched configura- 
tion is cooled 25 steps, the gluon content changes 
dramatically the fluctuations in the action 

and topological charge densities decrease two or- 
ders of magnitude, self-dual instantons and anti- 
instantons become apparent, ^ decreases from 
over 20,00 to 64, the string tension reduces to one 
fourth its original value, the lattice spacing de- 
fined by the nucleon mass decreases from 0.168fm 
to 0.142fm, and the magnetic hyperfine compo- 
nents of the quark-quark potential become essen- 
tially zero. Hence, for example, the energies and 



2.3. Zero mode zone expansion 

Because, as discussed in Section 3, cooling 
somewhat modifies the instanton distribution, it 
is useful to study the zero modes of the origi- 
nal uncooled configuration directly. Several stud- 
ies [||^0]-^2| calculate zero modes on the real axis 
and conjugate pairs displaced slightly off the axis 
due to the mixing of instantons and antiinston- 
tons, examine the discrete lattice counterpart of 
the index theorem, and consider the presently un- 
solved problem of separating physical continuum 
solutions from doublers. A clear indication of 
the role of zero modes in light hadron observ- 
ables is the degree to which truncation of the 

expansion ^^m+i^ ^° ^ e zer0 m0< ^ e zone 
reproduces the result using the complete propa- 
gator. The top plot in Fig. || shows the effect of 
including the lowest 16 through 128 of the 786,432 
modes on a 16 4 lattice in the vector correlation 
function, indicating that already 64 modes repro- 
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Figure 3. Truncated eigenmode expansion EC 



duce most of the strength of the rho peak and 
128 modes produce the full strength. The lower 
plot shows that 128 modes also reproduce the full 
pion strength |2(]] . Similarly, most of the strength 
of the disconnected graph contribution to the rf 
correlation function, which should be particularly 
sensitive to instantons, is already provided by the 
lowest 32 eigenmodes, and the fermionic defini- 
tion of the topological charge is nearly saturated 



by the lowest 8 modes |21 



2.4. Localization of zero modes 

The observation of localization of low eigen- 
modes of uncooled configurations at the locations 
of instantons identified by cooling confirms sev- 
eral important ideas. Even in the presence of fluc- 
tuations several orders of magnitude larger than 
the instanton fields themselves, the light quarks 
essentially average out these fluctuations and pro- 
duce localized peaks at the topological excitations 
as expected from semiclassical arguments. In ad- 
dition, it is clear that the instantons which sur- 



Figure 4. Isosurfaces on a periodic lattice 
of topological charge density FF(x) calculated 
from a cooled configuration (dark surfaces) and 
ipo(x)ipo(x) calculated on the corresponding un- 
cooled configuration (light surfaces) |B3[. 



vive cooling represent topological structure in the 
original configuration. This localization has been 
observed for unquenched Wilson fermions at T 
= 0, = 5.5 [^0| and for unquenched staggered 
fermions at t < T c at (3 — 5.65 ^3|. Isosur- 
faces of , ipo(x)ipo(x) for the lowest uncooled eigen- 
mode [£3| are compared with the cooled topolog- 
ical charge density in Fig. |4j 

2.5. Molecular correlations 

At high temperature, the fermion determi- 
nant favors alignment of instanton-anti-instanton 
pairs in the thermal direction pC| j. Although 
aligned molecules have not been observed directly 
in the topological charge density of cooled config- 
urations, a suggestive anisotropy in the topologi- 
cal susceptibility measured on differently oriented 
sublattices was reported at this conference |P3[ , 
Lattice slices that are narrow in the time di- 
rection but span the full spatial volume, which 
would have a large range in the value of Q for 
a dipole layer, were observed to have (Q 2 ) up to 
twice as large as slices that are narrow in one spa- 
tial direction but span the full time range, which 
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would tend to have canceling positive and nega- 
tive charges from the dipoles in the slice. Whereas 
it is not presently understood why a comparable 
asymmetry is also observed in quenched SU(2), 
subvolume susceptibility measurements are an in- 
teresting and promising development. 

2.6. Spectral density and (tpip) 

Calculations of the lowest 8 eigenvalues for a set 
of MILC configurations have been used to evalu- 
ate the density of states at low virtuality, p(0), 
and thereby numerically check the Banks-Casher 
relation on the lattice p3| . The results of low A 
fits to the integrated histograms yield values of 
(-0'0)m,-^o (to be compared with yp(0) in paren- 
theses) of 0.0114 (0.0138) at /3 = 5.65 and 0.002 
(0.002) at (3 = 5.725. 

The significant developments in the use of do- 
main wall chiral fermions have resulted in a num- 
ber of studies of zero modes reported at this 
conference [Q. For sufficiently large separa- 
tion between domain walls, there is an integer 
index for the lattice Dirac operator correspond- 
ing to the winding number for smooth fields as 
well as exact zero modes |2j|. In the presence 
of zero modes, the chiral condensate becomes 

-m = ^Ea =0 ^ + vL>oiw. leadin § 

to a ^ mass dependence. Recently [ p6| , calcula- 
tions with domain walls separated by 10 on a 16 4 
lattice in the presence of a smooth configuration 
with winding number one plus random fluctua- 
tions demonstrated this — dependence over the 

m 

range 10~ 5 < m < 10~ 3 . 

3. INSTANTON CONTENT OF THE 
QCD VACUUM 

Several recent studies of the instanton content 
of the SU(3) vacuum provide an unprecedented, 
albeit imperfect, view of the size distribution of 
instantons and their spatial correlations. To ap- 
preciate the issues involved in determining the in- 
stanton content, it is useful to briefly review the 
methods that have been used. 

3.1. Methods to extract the instanton 
content 

Standard cooling jl5],[l6) applied to QCD mini- 
mizes the action locally at each link. For the Wil- 



son action, one calculates Maxju] (Re Tr UJ] sta P les ) 
so that U _1 is the projection on SU(N) of 
S sta P ic S= ^ 3 , In the case of SU(2), since the 
sum of staples is a multiple of a group element, 
cooling makes the replacement j — > ^ M 3 ■ F° r 
SU(3), it is convenient to sequentially minimize 
with respect to SU(2) subgroups J27]]. 

Underrelaxed cooling [^8| takes a linear combi- 
nation of the original link and the sum of staples, 
\ — > (1 — c) I + I J2n 3 > which corresponds to 
APE smearing |29| and for small c has the ef- 
fect of making the minimization independent of 
the order of local updates. This has been used 
in Ref. |3(| with c =0.86 and sequential applica- 
tion to SU(2) subgroups to analyze the UKQCD 
SU(3) configurations. This work also calibrates 
the cooling at each mixing parameter c and f3 to 
define a number of calibrated sweeps which pro- 
duces approximately equivalent evolution of in- 
stantons and instanton-anti- instanton pairs. In- 
stantons are identified by peaks in the topological 
charge density, Q(x), assuming additivity and fil- 
ters are applied to assure the appropriate shape 
and to avoid counting fluctuations. 

Renormalization group cycling, as origi- 
nally implemented Ipjjj , is based on the 
renormalization group transformation e s ( u ^ = 
J d[V]e s ^ +KT ^ uy \ The blocking transforma- 
tion kT(U, V) maps a configuration V on a fine 
lattice to a configuration U on a lattice twice as 
course, and in the weak coupling limit the inverse 
transformation is given by the saddle point condi- 
tion Min {v} [S(V)+KT(U, V)]. Making use of the 
fact that V is the smoothest configuration that 
blocks to U and that the renormalization group 
preserves the instantons and large scale structure 
of the configuration, the cycling procedure first 
inverse blocks from the course lattice to the fine 
lattice and then reblocks back to a course sublat- 
tice shifted one lattice spacing along the diagonal. 

Although cycling has many conceptual and 
practical advantages, including stabilizing instan- 
tons with p > 0.94a and smoothing configura- 
tions sufficiently to identify instantons before the 
string tension has been severely diminished, it 
is very expensive computationally. Hence, effi- 
cient parameterizations of the cycling transfer- 
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mation were studied |3^] in SU(2). After inves- 
tigating 56 combinations of paths up to 9 links, 
it was determined that a good parameterization 
was obtained by APE smearing, providing an 
additional rationale for this procedure. In ad- 
dition, for reasons that are not yet clear, this 
work showed that although the sizes of smooth 
instanton configurations are stable under APE 
smearing, instantons and anti-instantons found 
in a configuration generated by the Wilson action 
may change size (generally growing) linearly with 
the number of steps. In Ref. p3| ], APE smearing 
I -> [(1 - c) \ + f E M 3 }su(3) with direct pro- 
jection onto SU(3) was used with c = 0.45 and 
15-30 smearing steps to study the vacuum distri- 
bution of instantons. For comparison, UKQCD 
used c = 0.86 and 23-46 steps and the calculations 
of topological charge by the Pisa group discussed 
below used c = 0.9 and 2 steps. 

Instead of sequentially minimizing the Wil- 
son action, improved cooling ]34| ] minimizes an 
improved action constructed from nxm plaque- 

ttes, S — X^{n,m} c n,m£>n,m,-> where Sn,r, 
1 



& — J2{n,m}' 

E w Tr(l-C 



m "| ) . Up to five terms are 
included, corresponding to lxl, 1x2, 1x3, 2x2, 
and 3x3 plaquettes, and the coefficients c n , m 
are determined to eliminate leading corrections 
in a and also to create a slight barrier to pre- 
vent shrinkage of instantons above a critical size. 
Whereas the Wilson action for a smooth instan- 
ton monotonically decreases with p, leading to the 
eventual shrinkage and disappearance of instan- 
tons after hundreds of cooling steps, the five-term 
improved action is extremely flat and preserves all 
instantons greater than p c =2.3a indefinitely. A 
comparably improved topological charge opera- 
tor yields nearly integer charge after 5-10 cooling 
steps. Improved cooling has also recently been 
used to study the SU(3) vacuum p5|| . 

A cooling method |?(| that is particularly effi- 
cient for a parallel computer is simultaneous re- 
laxation of the link variables at all lattice sites by 
discretizing the relaxation equation ^ = — 
Parameterizing an SU(3) matrix in the vicinity 
of U with by U[F] — e lF U, where H is hermitian 
and traceless, and taking a small step At yields 
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Figure 5. Sketch of extrapolations of the average 
instanton size to n c = cooling steps for data 
from Refs. @ (left) and @ (right). Solid and 
open symbols denote quenched and full QCD re- 
spectively. 

where Us denotes the staple (or generalized staple 
for an improved action) associated with U. Suffi- 
ciently small At produces cooling histories com- 
parable to overrelaxed cooling, and At = 0.025 
was used to analyze SU(3) configurations. 

3.2. Vacuum instanton distribution 

Key features and salient results of five studies 
of instantons in the SU(3) vacuum are summa- 
rized in Table |l|. 

A primary focus has been measuring the dis- 
tribution of instanton sizes, n(p), and the corre- 
sponding average size p. The fundamental limita- 
tion is the fact that the distribution is modified to 
some extent by each cooling procedure. Although 
in principle, cooling with the Wilson action even- 
tually allows instantons to shrink and eventually 
fall through the lattice, for the small number of 
cooling steps used with it in these studies, the 
effect is inconsequential. A serious problem for 
all methods, however, is the fact that sequential 
minimization of the action will rotate the rela- 
tive color orientation of instanton-anti-instanton 
pairs into the attractive direction and bring them 
together until they annihilate. Even over the rela- 
tively small range of cooling steps in Ref. , for 
example, instanton number densities range from 
55 fm~ 4 to 3 fm~ 4 , and Refs. 0j3^] extend down 
to 0.3 fm~ 4 . Continued evolution would even- 
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Table 1 



Studies of the instanton content of the SU(3) Vacuum 



a 


Lattice 


Method 


Pn/v 
(im) 


Pextrap 

(tmj 


N 
V 

(tm ) 


Ref. 


6.0 

6.2 

6.4 


16 3 x 48 
32 3 x 64 
24 3 x 48 
32 3 x 64 


Underrelaxed 
Cooling 


0.60(5) 


0.37(5) 


55-3.2 


ejj 


5.85 

6.0 

6.1 


12 4 
12 4 ,16 4 
16 4 


APE Smearing 




0.32 c 


1.1 


§§ 


5.85 
6.0 


12 4 
16 4 


Improved 
Cooling- 


< 0.53(5) tl 




3.3-0.38 




5.7 


16 3 x 24 


Cooling 




> 0.39 e 


0.59-0.28' 


m 


5.85 


16 4 


Relaxation 


0.50(5) 9 


0.43(5)^9 


5.3-1.49 




unquenchcd, k=0.16 
5.5 1 16 4 




0.52(5) 9 


0.42(5) fc ^ 


6.5-1.89 




Summary 


0.54(5) 


0.39(5) 





a Value 0.56(5) at N/V= 8.5 for = 6.4 evolved to N/V=3.2 using /3 = 6.0 data. 

b Extrapolation sketched in Fig. km 

c Value 0.3 scaled 5.6% using a(^fe = 440 MeV). 

d From graphs of N/V = 1.81 and 1.43 data. Evolution to 3.2 would reduce p further. 

e Value 0.36 scaled 9% using a(yfa = 440 MeV). Correlation function range underestimates average p |34| 

f Lattice spacing increased 9% using a(yfo = 440 MeV). 

9 Lattice spacing from hadron masses increased 18% using a(^/a = 440 MeV). 



tually annihilate all the pairs until only instan- 
tons or anti-instantons remain. The trade-off in 
all measurements is the uncertainty in identifying 
instantons and anti-instantons in the presence of 
fluctuations for a small number of steps versus the 
damage to the original instanton ensemble caused 
by a large number of steps. 

To compare the various calculations and em- 
phasize the similarity of their essential features 
as much as possible, two alternative definitions of 
the average instanton size have been used in Table 
|l|. One definition is the density extrapolated in 
the number of cooling steps, n c , to ji c = 0. This 
is more representative of the original, unfiltered 
vacuum but more subject to pattern recognition 
errors. As mentioned previously, it is necessary 
for APE smearing where instanton sizes were ob- 
served to evolve linearly. Data from Refs. |20 3(]] 
are plotted in Fig. |^ with linear extrapolations 
and the error bars bracketing these extrapolations 



are listed in Table 0. The calculation in Ref. fl7jl 
is effectively an extrapolation, since the scale cal- 
culated from hadron masses evolves nearly lin- 
early with the number of cooling steps. 

The second definition is the continuum limit of 
comparably filtered configurations. Ref. |5(J de- 
fines a calibrated number of cooling sweeps n c ((3) 
to reach the same density of instantons y at each 
(3 and extrapolates linearly in a 2 to the continuum 
limit . This is consistent with the fact that observ- 
ables such as the instanton distribution shown in 
Fig. U agree well at (3 = 6.0, 6.2, and 6.4. 

To facilitate comparison in Table [l], as indi- 
cated in the footnotes, authors' data were scaled 
to a lattice spacing determined by the string ten- 
sion Ja — 440 MeV and on was evaluated as 

close to y — 3.2 fm -4 as possible. From the table, 
there is qualitative agreement among all the cal- 
culations and the extrapolated value /3 ex trap li es 
in the range 0.39 ±.05 fm, consistent with the in- 
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Figure 6. UKQCD instanton size distribution 
30 1 . In top plot, p = 1 corresponds to 0.45 fin. 



stanton liquid model, while the value after cooling 



to -pr on the order of 3.2 m yields p« substan- 
tially larger, in the range 0.54 ±.05 fm. Although 
no attempt has been made to calculate w charac- 
terizing the uncooled vacuum, it is clear that the 
density of instantons identified during the early 
stages of cooling is substantially larger than the 
value 1 fm~ 4 in the instanton liquid model, which 
should likely be thought of as an effective density 
in the size range dominating physical processes. 

Each of the studies in Table |l| also calculated 
the distribution of instanton sizes, with similar 
results. The UKQCD results @, which had 
the highest statistics and greatest range of pa- 
rameters, are shown in Fig. ^. The upper plot 
shows the high level of consistency obtained be- 
tween comparably cooled configurations with (3 
= 6.0, 6.2, and 6.4 for 23, 46, and 80 cooling 
steps denoted by open circles, closed circles, and 



Figure 7. Distribution of like (solid circles) and 
opposite (open circles) charges as a function of 
distance R from a reference charge where 
lfm corresponds to R = 2.2 

crosses. Similar consistency was demonstrated at 
[3 = 6.0 between configurations on 16 3 x 48 and 
32 3 x 64 lattices. Full QCD with ^ = 0.82, has 
a similar distribution poj . The log-log plot |36| 
of the UKQCD data shows that the distribution 
of instantons grows roughly as p 6 at small p as 
expected in the dilute approximation, falls sig- 
nificantly below p e in the region of 0.3 fm, and 
decreases at large p much more rapidly than p 5 
characterizing a frozen coupling constant. 

Interesting correlations between instantons 
have also been observed |3(J] . Fig. |t] shows the dis- 
tribution of like and unlike topological charges as 
a function of the distance from a reference charge 
and clearly reflects the effect of instanton-anti- 
instanton attraction. Quantitatively, the average 
distances to the nearest like and unlike charges 
are 0.49 and 0.45 fm respectively, to be com- 
pared with the average size 0.56 fm. In a dilute 
gas, the distribution would be Poisson implying 
(Q 2 ) = Nj + Na , which only holds for p < 0.5fm. 
By all measures, the instanton distribution is not 
dilute, and is best thought of as a liquid. In 
addition, smaller and larger instantons tend to 
have opposite charge, with the smaller instantons 
having the sign of the total Q, so that the large 
charges are overscreened. Many other interesting 
details are found in fl30tl. 
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Table 2 



Topological Susceptibility 



X 1/4 (MeV) Method Rcf. 


SU(2) 


230 (30) RG Cycling 
220 (6) APE Smearing 
200 (15) Improved Cooling 
198 (8) APE + Renorm. 

Geometric + Renorm. 
226 (4) Spectral Flow 


31 

32 
34 

3b 
39 


1 
\ 
\ 


SU(3) 


187 (14) Underrelaxed Cooling 
192 (5) APE Smearing 
185 (9) Improved Cooling 
175 (5) APE + Renorm. 
197 (4) Spectral Flow 


3( 
33 
3C 
4( 
39 


] 
} 
} 
} 


180 Veneziano-Witten 


% 


§1 



3.3. Topological susceptibility 

It appears that we now understand how to de- 
fine quantities on a lattice that correspond to 
the topological charge in the continuum limit 
and that, as shown in Table [^, all methods for 
calculating the topological susceptibility are rea- 
sonably consistent with each other and with the 
Veneziano-Witten formula. 

As emphasized in Ref. |38|, all measure- 
ments of the susceptibility involve both addi- 
tive and multiplicative renormalization, a A \ l = 
^(/3)- 2 (xLttice- M(/3)). By placing known 
smooth configurations on a lattice and heating 
by Monte Carlo updates, M was measured for 
Liischer's geometrical charge and both M and Z 
were measured for the clover approximation to 
FF with zero and two APE smearing steps, lead- 
ing to a consistent value of x f° r a U three calcu- 
lations pq] . Cooling automatically removes the 
additive term and brings Z to unity by removing 
quantum fluctuations. It determines the suscep- 
tibility accurately because the total charge is un- 
affected by pair annihilation and because, with 
improved charge operators, the charge converges 
to a stable integer after very few cooling steps. 
Hence, the various cooling methods yield consis- 



tent results as shown. Finally, determination of 
the topological charge from spectral flow |39| also 
appears to be consistent with other methods. 

The susceptibility in SU(3) has also been stud- 
ied using the geometrical charge |37j . The Wilson 
action yields x 1 = 228 MeV, which is expected to 
be high because of the additive contribution from 
dislocations. Calculations with a renormalization 
group improved action that suppresses disloca- 
tions appear to have less additive renormaliza- 
tion, but lack of an accurate measurement of the 
scale precludes quantitative comparison. 

3.4. Equilibration of topological sectors 

An important problem in Monte Carlo calcula- 
tions with dynamical fermions is ergodically sam- 
pling the relevant topological sectors. Recent hy- 
brid Monte Carlo calculations @ with 3000 to 
5000 trajectories on 16 3 x32 lattices using Wilson 
fermions have demonstrated sufficient tunneling 
to achieve equilibration of Q at ^ 2L = 0.84, 0.76, 
and 0.69. At the lightest quark mass, an approx- 
imately Gaussian distribution of Q is observed 
and the integrated autocorrelation time for Q, r = 
54(4), is only moderately larger than the value 
r=42(4) for other observables. For ^ = .56, the 
mean value of Q did not equilibrate to within 
3,500 trajectories. These results are consistent 
with previous HMC calculations with staggered 
fermions, for which autocorrelation times of the 
order of 10 were observed for ^ =0.75 and 0.65, 
but ^ = 0.57 exhibited little mobility over 450 
trajectories f42]| . Thus, there is an intermediate 
range of masses down to — ~ 0.69 for which 
equilibration of topological sectors is achieved for 
both Wilson and staggered fermions. 

4. INSTANTONS AND THE STATIC 
POTENTIAL 

Several puzzles concerning the contribution of 
instantons to the static potential have emerged 
in recent years. A provocative calculation |4j| 
suggested that an instanton liquid with a distri- 
bution ~ p~ 5 at large p produced a linear con- 
fining potential with the physical string tension. 
In addition, a confining potential appears to sur- 
vive renormalization group cycling for SU(2), but 
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largely disappear when the instantons identified 
in the cycled configuration are replaced by in- 
stantons of comparable size and random orien- 
tation @. 

These and related issues are addressed by high 
statistics calculations [Q of the potentials in an 
instanton liquid with fixed instanton size and 
with the size distribution p 6 (/Og' 5 + p 3 - 5 )~~ , cor- 
responding to a frozen coupling constant . The 
potentials are linear at small distance and ap- 
proach a constant at large distance as expected 
analytically Ej|. With p = ^ fm and y = lfm~ 4 , 
both distributions yield the same slope, a = 0.11 
GeV/fm, roughly one tenth of the physical string 
tension. Whereas the fixed size potential ap- 
proaches a constant already by 1.5 fm, the p~ 5 
potential is still rising substantially at 3 fm due to 
the contributions of large instantons. Because the 
potential is roughly proportional to y and since 
by scaling, the physical string tension would be 
obtained at y = 1 fm -4 by increasing p by 10^ 
to 0.59 fm, one can understand several other re- 
sults. The peak of the distribution in Ref. (is) is 
at 0.4 fm corresponding to p above 0.5 fm, so the 
slope should be near the physical value. Also, the 
potential only extended to 1 fm, so the departure 
from confinement was not observed. Similarly, 
cooled distributions with larger instanton density 
and p are also enhanced. 

The question of what correlations remain in 
the renormalization group cycled configurations 
yielding cither confinement or at least a large 
slope out to 1 fm that are not contained in a 
comparable random distribution of instantons is 
interesting and should be pursued. One sugges- 
tive feature is the fact that an instanton liquid as 
usually constructed does not have Z(N) symme- 
try, whereas confining configurations generated 
by sampling the action automatically have it. 

5. INSTANTONS AT NONZERO 
TEMPERATURE 

5.1. Topological susceptibility 

Measurements of the topological susceptibil- 
ity beyond the critical temperature confirm the 
idea that instantons should be suppressed by De- 
bye screening, where the leading suppression fac- 



tor {H is e -^ {2N " +N ' ){ ' npT) ' 2 As at zero tempera- 
ture, within one's ability to read numbers from 
graphs, quenched measurements using one and 
two APE smearing steps corrected for additive 
and multiplicative renormalizatio n [ffO| agree with 
analysis of cooled configurations |47|. The cooled 
measurements showed that the average instanton 
size decreased from p = 0.33 fm at low temper- 
ature to 0.26 fm at 1.34 Tc and that with the 
latter value, Debye screening fit the rapid fall off 
of the susceptibility. Comparable cooling calcu- 
lations with staggered fermions |Q showed anal- 
ogous results, with the instanton size p = 0.44 fm 
at .75 Tc decreasing to 0.33 fm at 1.3 Tc and 
Debye screening with the latter size fitting the 
susceptibility above Tc- The transition is con- 
siderably sharper than in the unquenched case, 
consistent with the N / dependence of the Debye 
screening. 

5.2. Calorons 

Gauge fields at nonzero temperature may 
be classified by the Polyakov loop Voo = 
LinXc^ooPe-'o dtA o( x ' t ) ^ anc j one usua lly considers 
Voo = ±1 since fluctuations around classical so- 
lutions with other values produce a non-zero vac- 
uum energy density Ji9"| . The periodic instanton, 
or caloron is given by A fJ _(x) — ^ff^ y T^d v ln(j), with 
(f> a periodic sum of instantons displaced in the 
thermal direction p0| . In the limit where the in- 
stanton size p is large compared with /?, the spa- 
tial distribution becomes independent of t, and 
gauge transformation yields a BPS monopole at 
the center of the instanton (and a monopole with 
opposite charge at infinity |plf . 

During the past year, a remarkable set of solu- 
tions comprised of monopoles has been discov- 
ered |32|,[53| for nontrivial values of Voo , cor- 
responding to to a periodic array of instantons 
twisted by Voo- The closed form solutions are 
surprisingly simple. For SU(2) withPoo = e 2 ™uir 3 
the action density is tcF 2 v = d 2 d 2 In ip, where ip = 

-cos(2^+cic 2 + r?+ g i + r ^ p % lS2 +V(^f+^f), 
the positions of two constituent monopoles are de- 
noted r m = \x-y m \ for m=l, 2, \y 2 - y\\ = np 2 , 
v x = 2uj, v 2 = 1 - 2lo, Cm = cosh(27r^ TO r m ), 
s rn = sinh(27w m r m ), and the solution is scaled 



12 




Figure 8. Action density profiles for calorons |5: 
at ui = 0, g, and \ with p = j3 = 1. 



to /? = 1. General expressions are given for tri 7 ^ 
and A^x) for SU(N) in Ref. f|. 

In the SU(2) case, for p large relative to (3 the 
the solution separates into two lumps separated 

2 o 2 

by with monopole charges ±1, masses sZL^m | 
and widths ~ In the limit of large 4 the 

lumps become well separated, spherically sym- 
metric, and time independent. Since they are self- 
dual, they are therefore BPS monopoles. As lo ap- 
proaches or \ , the second lump becomes negligi- 
ble and the solution becomes the spherically sym- 
metric Harrington- Shepard caloron. The action 
densities plotted on logarithmic scales in Fig || at 
t=0 show the separation into lumps at p = fl. For 
small ^ , the caloron approaches a single instanton 
solution. In the maximal abelian gauge, as p in- 
creases from to oo, the small monopole loops as- 
sociated with instantons grow. When these loops 
extend from to (3, they merge with their peri- 
odic extensions to produce pairs of monopole lines 
which eventually become the well-separated BPS 
monopoles |5^] . In the case of SU(N), the solution 
generalizes to N monopole lumps, and results are 
shown in Ref. |p2| . 

These novel finite temperature solutions pro- 
vide interesting examples of gauge fields with 
topological charge which are built out of 
monopoles without gauge fixing, and we may 
hope that they will lead us a step closer to ana- 
lytical insight into the role of monopoles in QCD. 



6. SUMMARY AND OPEN PROBLEMS 

We have come a long way in exploring instan- 
ton physics on the lattice. The behavior of two- 
point functions, the similarity of cooled and un- 
cooled calculations, and the contributions of zero 
modes provide substantial evidence that instan- 
tons play a major role in the physics of light 
hadrons. There has been significant progress in 
calculating quark zero modes, observing local- 
ization, extracting the instanton content of the 
vacuum, calculating the topological susceptibility, 
and discovering analytic solutions with monopole 
constituents of instantons. Salient problems for 
further investigation include reducing the ambi- 
guity in the instanton distribution in the vacuum, 
removing the effect of doublers in spectral flow 
analyses, understanding the role of instantons in 
confinement, and especially studying full QCD 
with light chiral fermions, where the contribu- 
tions of instantons and zero modes become most 
significant. 
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